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The un-integrated parton distribution functions (uPDFs) obtained from a CCFM 
evolution are studied in terms of the intrinsic transverse momentum distribution 
at low scales. The uPDFs are studied for variations of the renormalization and 
factorization scales. 



1 Introduction 

Un-integrated parton distributions are best suited to study the details of the per- 
turbative as well as the non-perturbative QCD evolution. The general form of the 
integral equation for the parton evolution is: 

f dz f dq 2 ~ / x \ 

xA(x,k t ,q) = xA (x,k t ,q) + J — J <d (order) A s P(z, q, k t )xA y-, k' t , qj (1) 

with q being the evolution (factorization) scale. The ordering condition of the 
evolution is specified by 8 (order) which in the case of CCFM is given by 8 (order) = 
Q(q — zq). The first term of the rhs of eq.Q gives the contribution of non-resolvable 
branchings between the starting scale Qq and the factorization scale q and is given 
by: 

xAo(x,k t ,q) = xA (x 1 k t ,Q )A s (q,Qo) (2) 

where the Sudakov form factor A s (q2, qi) describes the probability of no radiation 
between Qq and q. The second term of the rhs of eq.Q describes the details of the 
QCD evolution, expressed by the convolution of the splitting function P with the 
parton density and the Sudakov form factor A s . 



2 The starting distribution xAq(x, kt, Qo) 

The non-resolvable branching contribution is directly related to the initial (starting) 
distribution Aq(x, kt , Qo), and especially to the intrinsic kt distribution. Often a 
Gaussian type distribution with width fc is assumed and Ao can be parameterized 
as: 

xA {x,k t ,Q ) =Nx p °(l-x) Pl ■ exp (-kf/k 2 ) (3) 

where N is a normalization constant, and po,Pi being parameters to be determined 
experimentally. 

In the following sections the dependence of the choice of po,Pi and fco will be 
investigated. The CCFM evolution equation is applied to the starting distribution 
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xAq(x, k t , Qo) and then convoluted with the off-shell matrix elements. The simple 
gluon splitting function is applied |1I2| . 

„ , , % <5s(<7 2 (l - z) 2 ) a s (k?) , , , . . A . 

P gg {z,q,k t ) = {_ z + ^^A ns (z,q 2 ,k t ) (4) 

with q = p t /(l — z) and the non-Sudakov form factor A ns is defined by: 

logA ns (z i ,<?,k t ) = -a s J ^ y ^ • 9(fc t - g)0(g - z'q) (5) 

The collinear cutoff Q g — 1.3 GeV regulates the region of z — > 1, applied both to 
the real emissions as well as in the Sudakov form factor for the virtual corrections, 
which is equivalent to using the plus prescription of the splitting function P + . 

Due to the angular ordering in CCFM a kind of random walk in the propagator 
gluon kt can be performed, and values of small kt can be reached. A cut k^ ut is 
introduced to avoid this region. 

The CCFM evolution equations have been solved numerically PQ using a Monte 
Carlo method. In Tab. ^ the fit results for the un-integrated gluon densities are 
summarized. The input parameters are varied such that after convolution with the 
off-shell matrix elements the results agree best with the measured structure function 
F 2 (x, Q 2 ) in the range x < 5 • 10" 3 and Q 2 > 4.5 GeV 2 as measured at HI [311) 
and ZEUS |5lf>| . The starting distribution is given by eq.®, where the parameter 
pi = 4 is kept fixed, but po and the width of the Gaussian ko was adjusted to give 
the smallest x 2 /points. The results are given in Tab-^and labeled as set AO. 

Equally good fits can be obtained using different values for the soft cut k£ ut and 
a different value for the width of the intrinsic k t distribution set BO. Motivated 
by the success of the saturation model |7], the parameterization of GBW can also 
serve as a starting distribution: 

xA (x, h,Qo) = Ro( x ) k t ex P (- R o( x ) k t) ( 6 ) 

with the parameters given in 7 . As seen from Tab. ^a reasonable description of 
the data can be achieved set GBW. 

The renormalization scale dependence of the final cross section can be estimated 
by changing the scale used in a s in the off-shell matrix element. Since here we are 
using the LO a s matrix elements, any scale variation will change the cross section. 
In order to obtain a reasonable result, the uPDF was fitted to describe F 2 by 
varying the scale /i r . The set A- (set B-) correspond to a scale /i r = 0.5p t and 
the set A+ (set B+) correspond to a scale \l t = 2p t . by: 0.5 • /i r and 2 • \i r . 
The resulting x 2 /points is also shown in Tab. ^ However changing the scale \i r 
and /i/ simultaneously with the same factor results in a much larger % 2 , which is 
understandable, since the length of the evolution ladder is changed by [if. 

The CCFM evolution is performed in an angular ordered region of phase space 
and the factorization scale is related to the maximum allowed angle for any emission. 
In CCFM the scale Uf (coming from the maximum angle) can be written as: 



u f = A/* + Ql (7) 



Un-integrated Parton Density Functions in CCFM 



3 



with s being the invariant mass of the qq subsystem, and Q± its transverse mo- 
mentum. In this definition, the scale fif is related to the quark pair. 

In order to investigate the uncertainties coming from the specific choice of the 
evolution scale given in eq.(|7J), another definition is applied, relating the factoriza- 
tion scale only to the quark (or anti-quark) : (if = with p t being the transverse 

momentum of the quark(anti-quark) and z = ■ This definition follows closely 
the definition of the rescaled transverse momentum q in CCFM. In Tab.^is shown, 
that also with this definition a reasonable description of the data can be achieved, 
but the parameter po of the starting distribution is changed significantly. 



set 


M/ 


fir 


Po 


k 


k cut 


X 2 /points 






= \/Pt + m2 




(GeV) 


(GeV) 


X 2 /points 


AO 


+ Ql 


1 ■ H 





1.33 


1.33 


276/248 = 1.1 


A0- 


\/s + Ql 


0.5 • fif 


-0.01 


1.33 


1.33 


433/248 = 1.75 


A0+ 


\/s + Ql 


2-M/ 


-0.01 


1.33 


1.33 


316/248 = 1.3 


Al 


pt 




-0.1 


0.8 


1.33 


275/248 = 1.1 


BO 


V s + Ql 


1 ■ fif 





0.8 


0.25 


356/248 = 1.4 


B0- 




0.5 -fif 


0.01 


0.8 


0.25 


483/248 = 1.95 


B0+ 


V s + Q'i 


2- ii f 


0.01 


0.8 


0.25 


297/248 = 1.2 


Bl 


pt 

1-z 


1 


-0.1 


0.8 


0.25 


702/248 = 2.8 


GBW 


V s + Ql 


1 






0.25 


349/248 = 1.4 



Table 1. The different settings of the CCFM un-integrated gluon densities. The initial distribution 
is parameterized as xAq(x, kt, Qq) = Nx P( >(l — x) pi ■ exp (— k^/k^J with pi = 4 and Q g = 
1.33 GeV. In the last line the starting distribution from the parameterization of GBW is used. 
In the last column, the \ 2 /points to HERA F2 data I3I4I5I6I is given (for a; < 5-1CT 3 and Q 2 > 4.5 
GeV 2 ). 



3 k t dependence of xA.(x, kt, fif) 

In Fig.^the uPDFs (as described in Tab.^) are shown as a function of kt for differ- 
ent values of x. Also shown is the contribution of the first term of the rhs of eq. (|T|l . 
It can be seen, that the influence of the intrinsic k t distribution is concentrated 
at small values of kt, whereas at values kt ~ 1 GeV the perturbative evolution is 
responsible. The specific choice of the fef 4 ' results in different distributions at small 
values of k t . It is also interesting to observe, that the dip in k t visible in Fig. ^ 
(upper plot) is directly related to fe™*. 

It is interesting to observe, that even with very different starting distributions, 
the uPDFs after perturbative evolution are similar (at larger kt). Thus the small 
kt region provides new information on the non-perturbative part of the parton 
density functions, which can only be investigated by using uPDFs. Especially the 
question of saturation effects at small k t can be directly investigated with uPDFs. 
The difference of using a Gaussian intrinsic k t distribution to a form motived by 
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the saturation model of GBW is significant. It is also important to note, that 
these differences are washed out for inclusive distributions, as all uPDFs are able 
to describe F2 at a similar level. 

q = 10 GeV 




10 2 10 1 1 10 10 2 io~ 2 io' 1 1 10 10 2 io~ 2 io' 1 1 10 10 2 10 3 

*2 (GeV 2 ) 

Figure 1. Comparison of the different sets of un- integrated gluon densities obtained from the 
CCFM evolution as described in the text. 
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